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■ Abstract 

. We compute finite-size corrections to nucleon masses and magnetic moments in a periodic, spatial 

psj \ box of size L, both in QCD and in partially-quenched QCD. 
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I. INTRODUCTION 



Impressive progress is being achieved in deriving properties and interactions of hadrons us- 
ing lattice QCD. In several instances, lattice methods are making predictions of hadronic 
quantities at the several-percent level Despite remarkable technical advances, current 
computational limitations continue to necessitate the use of quark masses, rUq, that are sig- 
nificantly larger than the physical values, lattice spacings, a, that are not significantly smaller 
than the physical scales of interest, and lattice sizes, L, that are not significantly larger than 
the pion Compton wavelength 2]. Therefore, lattice QCD simulations of hadronic physics 
require extrapolations in the quark masses, lattice spacing and lattice size, and ultimately 
it is confidence in these extrapolations that will allow a confrontation between lattice QCD 
and experiment. Fortunately, in many cases, the dependence of hadronic physics on these 
parameters can be calculated analytically in the low-energy effective field theory (EFT). 
Calculability requires maintaining the hierarchy of mass scales, 

|p I , < < fl-i , (1) 

where |p | is a typical momentum in the system of interest, is the pion mass and 
~ 2v^vr/ is the scale of chiral symmetry breaking (/ = 132 MeV is the pion decay 
constant). In a spatial box of size L, momenta are quantized such that p = 27rn/L with 

G Z. The hierarchy of eq. then requires maintenance of the additional inequality 
f'L ^ 1. This bound ensures that (non-pionic) hadronic physics lives inside the box. In 
addition, the bound (m7rL)^(/L)^ ^ 1 ensures that the box size has no effect on spontaneous 
chiral symmetry breaking 0,0]. These two bounds, taken together, then imply that we must 
have m.„L> 1. When {mT^LY^f 1, and therefore m^L -C 1, momentum zero-modes 
must be treated nonperturbatively |^ and one is in the so-called e-regime. 

Here we will consider the range of pion masses, 139 MeV < m^r < 300 MeV, and therefore 
we will take L> 2 fm, keeping in mind that the EFT may be reaching the limits of its validity 
when this bound on L is saturated, particularly when the pions are light. For the observables 
considered here, finite-volume effects tend to be small for L > 4 fm. It is therefore of 
interest to have control over the finite-size dependence of hadronic observables in the range 
2 fm < L < 4 fm. Chiral perturbation theory (^PT), which provides a systematic description 
of low-energy QCD near the chiral limit, is the appropriate EFT to exploit the hierarchy of 
eq. (H)) and to describe the dependence of hadronic observables on 1, ^ ,1 ^ . Recent work 
has investigated the finite- volume dependence in the meson [3,0,E3,EI, 



— , TiI3HHHE3 

sector and in the baryon [18|, |l9|, |20|, |2l|, |2^ sector. In this paper we compute the leading 



finite-volume dependence of the nucleon masses and magnetic moments in baryon xPT, 
including the A as an explicit degree of freedom ^. The finite-size dependence of the nucleon 
mass was first studied in Ref. 0], and has recently been computed to 0{m\) in baryon xPT 
(without including the A as an explicit degree of freedom) in Ref. j2l[ . (Some discussion of 
the effects of the A on the finite-size dependence of the nucleon mass appears in Ref. jl8l|.) 

We also give expressions for the finite-size dependence of the nucleon masses and magnetic 
moments in partially-quenched QCD (PQQCD), including strong isospin breaking. The cost 
of simulating dynamical quarks with light masses suggests separately varying the sea and 



Recent work has suggested that for certain observables, a rearrangement of the chiral expansion 

may improve convergence. We do not utilize these modified chiral expansions in this paper. 
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valence quark masses in the lattice QCD partition function, a procedure known as partial- 
quenching. This procedure has important advantages beyond issues of cost; by increasing 
the dimensionality of the parameter space that is explored, lattice QCD simulations can 
provide additional "data", which can significantly improve the qualit y o f extrapolations. 
XPT has been extended to describe both quenched QCD (QQCDjJlpJipJp, H Q with 
quenched chiral perturbation theory (QxPT) and PQQCD IsO. ElL l32l 33 with partially- 
quenched chiral perturbation theory (PQxPT). Recently, meson and baryon properties have 



been studied extensively in both QxPT M] and PQxPT [34, M |37|, lag. The 
effective field theory (EFT) describing the low-energy dynamics of two-nucleon systems and 
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nucleon-hyperon systems in PQQCD has also been explored 

This paper is organized as follows. In Section |nj the leading finite-size corrections to 
the nucleon masses are computed. The same is done for the nucleon magnetic moments in 
Section lllll We conclude in Section |TV1 Mathematical details and the partially-quenched 
extensions of the QCD results (including strong isospin breaking) are left to Appendices. 



II. THE NUCLEON MASSES 

A. The Infinite- Volume Limit 

For purposes of setting notation, we will begin by reviewing the derivation of leading terms 
in the chiral expansion of the nucleon mass. The relevant leading Baryon mass operators in 
two-flavor xPT are 

L = iNvVN + 2aMNM+N + 2aMNN tr [M+] - iT% -VT^ + A T^T^ , (2) 



where the chirally-invariant mass operator is M.^ = | + C'^gC); with 

diag(mu,md), and ^ = exp {iiiaTa/ f) is the usual two-fiavor Goldstone matrix. The rele- 
vant leading axial operators are 

C = 2gANS^A^N + g^^ [T'''" A^^, N, + h.c. ] . (3) 

The mass of the i-th nucleon has a chiral expansion of the form 

M, = Moifx) - Ml'\ii) - Mp)(/i) + ... , (4) 

where a term m/"'* denotes a contribution of order m^, and i = p,n. The nucleon mass is 
dominated by a term in the xPT Lagrangian, Mq, that is independent of m„. Here A, the 



A-nucleon mass splitting, is assumed to be of the same chiral order as the pion mass |42l |4 



Each of the contributions in eq. Q depends upon the scale chosen to renormalize the theory. 
While to O (mq), the objects Mq and M^^^ are scale independent, at one-loop level, 0(mq/^), 
they are scale dependent. The leading dependence upon niq, occurring at O {rriq), is due to 
the operators in eq. (0) with coefficients a^/ and ctm- The leading non-analytic dependence 
upon ruq arises from the one-loop diagrams shown in Fig. ^ 

In isospin-symmetric QCD ^, with mu,md — m, one finds the nucleon mass at one-loop 



^ The nucleon masses, including strong isospin breaking, may be obtained by taking the QCD limit of the 
partially-quenched expressions given in Ref. [s^. 



3 




FIG. 1: One-loop graphs that give contributions of the form ~ ruq to the masses of the proton 
and neutron. A solid, thick-solid and dashed line denote a nucleon, A-resonance, and a meson, 
respectively. The solid-squares denote an axial coupling. 



order in the chiral expansion 4J|, 

Mn = Moifi) - 2m {aM + ^au) (/u) 
where 

F(m,A,/i) = (m^-A^) ( VA^ -m2 log 



3 2 3 I ^9\n 



'A - VA2 -m? + ie 



A + VA2 - m2 + 



Alog 



1 . 2 , / 

-Am log — 

2 \ /i^ 



(5) 



(6) 



and A is the A-Nucleon mass splitting. Here we have used dimensional regularization (dim 
reg) with MS to define the divergent loop integrals and Mo(/i) and M'^^)(/i). 



B. Finite-Size Corrections 



In the infinite- volume limit, the nucleon mass may be written as 



.99', 



= Mo(7^)-2m(aM + 2aM)(7^) ~t^Moo,0) -z^J^(oo,A), (7) 



2/2 



where 



T7^(oo, A) 



3 Jn (27r)4 (A;o - A - te){ki - P - ml + ie) ' 
Here TZ denotes a choice of ultraviolet regulator and a renormalization scheme ^. 



In dim rcg with MS, 



'247r2 



F(m^,A,/i) , 



where F(to^, 0, fi) = 7rm^, which recovers the results of eq. 



(8) 
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In a spatial box of size L, X-ji generalizes to 

where we have rotated the integral to Euclidean space and accounted for the quantization 
of the momentum levels due to the periodic boundary conditions. Feynman parameterizing 
and explicitly evaluating the k/i integration lead to 

J7^(L,A) = -t- dX { ^ y ^—^ , (10) 

where = + 2AA + m^. We can now write the finite-size corrections to In as 
6lI{A) = In{L,A) - X7^(oo,A) 



6 Jo 



dX 



(11) 



^'\[k'+/3i]y') ^^^'[[k'+Pir/\ 

where 6l (/(I^I)) is defined in eq. flAl|) of the appendix. Notice that 5iX(A) is a purely 
infrared quantity and, as such, is independent of 71^. Using eqs. ((Tj) and (fTT|). the finite-size 
corrections to the nucleon mass may then be expressed as 



6lMn = Mm{L) - Mm{oo) = ^ ^lI{0) - i 5lJ(A) . (12) 

Using the "master" formula, eq. ()A6|) . derived in the appendix, we find 



5aiA) = -^/C(A), (13) 



where 

/C(A) = / rfA/?A ^ \{L\n\)-'K^{Lf3^\n\) - /?a i^ol^^/^A |n|)) . (14) 

Here Kn{z) is a modified Bessel function of the second kind. With A = the integral over 
A can be carried out explicitly (see Appendix) and one has 

/C(0) = -^miy{L\n\)-'exp{-L\n\m^). (15) 

Notice that this function contains no power-law corrections. Finally, we have 

^9a _ 9an 

This is the exact formula for the finite-size corrections to the nucleon mass at leading order 
in baryon xPT. In Fig. |2] (right panel), the ratio of the nucleon mass size dependence to 
the (infinite volume) nucleon mass has been plotted against the box size L for various 
pion masses. The solid and dashed lines correspond to the QCD formula of eq. (fTBj) with 
m^r = 139 MeV and 300 MeV, respectively. The dotted hne corresponds to PQQCD in the 
isospin limit taken from eq. ()B7|) with m^r = 139 MeV and = 300 MeV. We use the 
parameter set: / = 132 MeV, qa = 1-26 and g^^ = 1.4. 



= m - ^(A) • (16) 



^ This implies that finite-volume effects should be independent of the lattice spacing a, which appears 
implicitly as the ultraviolet cutoff 7r/a in all sums and integrals. 
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FIG. 2: Left panel: The ratio of the asymptotic formula, eq. jlTj ), over the exact formula, eq. 
as a function of L for various values of m-,^. The solid, dashed and dotted lines correspond to 
rriTT = 100 MeV, 200 MeV and 300 MeV, respectively. Right panel: The ratio of the nucleon mass 
size dependence to the (infinite volume) nucleon mass vs. L. The solid and dashed lines correspond 
to QCD with = 139 MeV and 300 MeV, respectively. The dotted line corresponds to PQQCD 
with = 139 MeV and mi = 300 MeV. 



C. The Asymptotic Limit 



Using eqs. ()15|1 and (jAllll . in the large-L expansion one has 

/ 9g\ml , Agl^rml/'^ 1 \ 1 



+ 



(27r)3/2/2 ALV2 j L 



exp {—mj^L) 



(17) 



where (5iMAr— 5^ Mat = (9(exp (— m.n-L)/L^/^). In the M^r — * cxo limit, the leading term in the 
large-L expansion is in agreement with Ref. [i^ and in disagreement with Ref. \ Fig.H 
(left panel) plots the ratio 6^Mn/SlMn- as a function of L for various pion masses. Clearly 
the utility of eq. (fT7|) is purely aesthetic; even with heavy pions, the asymptotic formula, 
eq. ()17|) . is not accurate for L < 10 fm. This points to the importance of exponential 
corrections; for itIt^V^ 1 convergence of the momentum sums requires keeping terms with 
\n\ > 1, i.e. one must include corrections of 0{exp (— |n|m^L)). 



III. THE NUCLEON MAGNETIC MOMENTS 



The Infinite- Volume Limit 



With the finite-size corrections for the masses, it is straightforward to get the magnetic 
moments. The leading operators contributing to the nucleon magnetic moments are 

e 



4M; 



N 



(18) 



For a detailed discussion of this disagreement, see Ref. 
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where F^^ is the electromagnetic field-strength tensor, Mn is the physical value of the 
nucleon mass, is the isoscalar nucleon magnetic moment, /ii is the isovector nucleon 
magnetic moment and t^_^_ = | ( + C''""^^) • 

In isospin-symmetric QCD one finds the nucleon magnetic moment matrix at one-loop 
order in the chiral expansion 4^ 4^, i^l 



N 



(a) 



N N 



(b) 



N 



FIG. 3: One-loop graphs that contribute to the proton and neutron magnetic moments. A solid, 
thick-solid and dashed line denote a nucleon, A-resonance, and a meson, respectively. The solid- 
squares denote an axial coupling and the solid-circles denote a leading- order electromagnetic inter- 
action. 



fi = /io + /il Ts 



N 



4np 



2 I 2 77 

9 A ^7r+ + g 9 AN TT 



(19) 



The scale dependence is left implicit. The proton and the neutron magnetic moments are 
the diagonal elements of /i. The first term within the brackets is from Fig. ^a) while the 
second term is from Fig. ^b). The function Ti = J^{mi, A, fi) is 



7r^(m, A,/i) = V — log 



'A - VA2 - m2 + ie" 



^A + VA2 -m"^ + ie J 
Here again we have used dim reg with MS. In the limit A — > 0, JF(m, 0, fi) = m. 



(20) 



B. Finite-Size Corrections 

In the infinite- volume limit, the nucleon magnetic moments may be written as 

2 



fl = /io + /^l T-3 



9 A Jn{oo, 0) + - 9 AN Jn{oo, A) 



T3 



where 



J7e(oo,A) 



d 



dml 



Jr,(oo, A) 



Therefore, the finite-size corrections to fi are 



hp' = p{L) - /i(oo) 



9\ hJ{Q) + - 91n hJ{A) 



(21) 



(22) 



(23) 



Using eqs. (fT^ . (fT^ and (j^ . and the properties of modified Bessel functions, it is straight- 
forward to find 



247r2 



(24) 



where 



POO 

3^(A) = / dXY. [?,Ko{L(3^\n\) - {L(3^\n\) K,{L(3^\n\))] . 



With A = one has 
3^(0) = 

Finally one arrives at 



TT 

-- 



2{L \n\ 



exp( 



\n\ 



(25) 



(26) 



M 



N 



9\ 3^(0) + g 3^(A) 



^3 • 



(27) 



This is the exact formula for the finite-size corrections to the nucleon magnetic moments 
at leading order in baryon xPT. In Fig. H] (right panel), the ratio of the proton magnetic 
moment size dependence to the (infinite volume) magnetic moment has been plotted against 
the box size L. The solid and dashed lines correspond to the QCD formula of eq. (j27j) with 
rriT, = 139 MeV and 300 MeV, respectively. The dotted line corresponds to PQQCD in the 



isospin limit taken from eq. (jBlljl with 



139 MeV and mi = 300 MeV. 
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FIG. 4: Left panel: The ratio of the asymptotic formula, eq. \2^) . over the exact formula, eq. \21\j , 
as a function of L for various values of m-,^. The solid, dashed and dotted lines correspond to 
m-j^ = 100 MeV, 200 MeV and 300 MeV, respectively. Right panel: The ratio of the nucleon 
magnetic moment size dependence to the (infinite volume) nucleon nucleon magnetic moment vs. 
L. The solid and dashed lines correspond to QCD with m^^ = 139 MeV and 300 MeV, respectively. 
The dotted line corresponds to PQQCD with = 139 MeV and m% = 300 MeV. 



8 



C. The Asymptotic Limit 

Again using eqs. (fTS|) . (fT^ and (j^ . together with eqs. and ()A11|) . in the large-L 
expansion one has 

where the dots denote contributions of (9(exp (— 777,7^^)/^^/^). Fig. IH (left panel) plots the 
ratio 6f^fi/6Lfi as a function of L for various pion masses. The curves are similar to those of 
the nucleon mass in Fig. |2l as is the conclusion about the practical utility of eq. (|28|) . 



IV. DISCUSSION AND CONCLUSION 

It is hoped that in the near future lattice (PQ)QCD simulations of baryon properties will 
encounter the chiral regime, where the quark masses are sufficiently small to allow a mean- 
ingful chiral expansion in quark masses, box size and lattice spacing. It is likely that this 
regime has been encountered in recent work on heavy-meson s yste ms m. 

The results of this paper, together with the results of Refs. [35| and |48| , give the depen- 
dence of the nucleon masses and magnetic moments on the sea and valence quark masses 
and on the lattice spacing, a, and size, L, to leading order in the chiral expansion. We 
eagerly await lattice (PQ)QCD simulations within the chiral regime where this parameter 
space may be fruitfully explored. 
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APPENDIX A: FINITE SIZE CORRECTIONS 



1. The Master Formula 

We wish to evaluate 

As this difference is ultraviolet finite, we omit the label denoting the scheme dependence 
of the individual sum and integral. This expression has been evaluated in many places ^. 



^ References that have been of use to the author include Refs. lid l49l Is 
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Using the identity 



D- 



1 POO 

m) Jo 



r(m) 



one finds 



[f2 + _A/f2]m J (47r)3/2r(m) Jo 



drjT] 



m-5/2 —riM'^ 



3/2 



L3 



I 



Expressing the momentum as / = 2TTn/L and using the Jacobi identity |49|, 

,2 



leads to 



S{z) 



1 



E 



S(z) 



-sC-) 

z z 



"Tim) ttIjo 



V [P + A<2]m ; (47r)3/2r(m) ^Jo 

Performing the integral over r] then gives the "master" formula 
^ 1 \ 2"5-"'A^3-2m 



21m 



7^3/2 Y[m] 



(A2) 



(A3) 



(A4) 



(A5) 



E (LA^|n|)-^+-ir|_(L>l|r?|), (A6) 



where -ft'n(-2) is a modified Bessel function of the second kind. 

A well-chosen change of integration variable and the properties of modified Bessel func- 
tions allow one to write eq. (fT^ as 

/c(A) = E(^Nr 7^ rdie{e-mi+/s.r"'K,{Lm))- m 

We find no useful simplification of this formula in the general case. With A = one directly 
finds 



- ml J2 {L\n\m^)'2Ki{L\n\m^) 

^ n^O 



(A8) 



which gives eq. (fT3jl. 



2. The Asymptotic Limit 

In the large-L limit, using the expansion of the modified Bessel function for large argument, 
one finds from eq. ()A7|) . 

/C( A) = 3v^ ^ ^ {l^" 12 ie -ml + A^)-^^ exp(-L ) + • • • (A9) 
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where the dots denote contributions of 0{exp [—mT^L)/L^/'^). Observe that one can expand 
the integrand in powers of = — m^, 



= m^/^ ^ exp (-L m^) + ... (AlO) 

where the dots denote contributions of C(exp (— m^L)/L^). (Similar technology has been 
developed in Ref. in the context of the heavy- meson systems.) Plugging this back into 
eq. ()A9|1 one finds, in the asymptotic limit, 

/C(A) = -3v^m^/2 _i_exp(-Lm^) + ... . (All) 

where the dots denote contributions of 0{exp {—mT,L)/L^^'^). 



APPENDIX B: PARTIALLY-QUENCHED QCD 
1. Nucleon Masses 

We work in PQQCD including isospin breaking, but with electromagnetism turned off ^. The 
Lagrangian describing the interactions of Bijk (containing the nucleon) and Tijk (containing 
A), which transform in the 70 and 44 of SU{A\2)v, respectivelv, with the pseudo-Goldstone 
bosons at leading-order in the chiral expansion is j28| 



C = 2a (BS^'BA^) + 2(5 (BS^'A^b) + a/|c[ (t" A,b) + {BA.T'') ] . (Bl) 



Here the axial-vector field A^ is a six-by-six matrix. Matching to the QCD effective La- 
grangian of eq. © and to the additional operator 

C = giNS^N tr [ A J , (B2) 

one finds that at tree-level, 

4 1 2 1 

« = + -^1 , P = -gi -^9A , C = -QAN ■ (B3) 

The finite-size corrections to the proton mass are given by 

5lMp = [ M[/c(m„„0)+/C(m„rf,0) + 2/C(m,„,0) + 2/C(mz„,0) + 3e?^„,,„(0) ] 

onj \ on 



^ For details we refer the reader to Ref. 13 
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+ [ /C(m„„, 0) - 5/C(m„rf, 0) + 3/C(mjd, 0) + 2/C(mj„, 0) + 3IC{mid, 0) + 2}C{miu, 0) 

DVT 

+3^,.,,.(0) + 66?,„,,,(0) + 36?,,,,,(0)] 

+ [ /C(m,„, 0) + /C(mi„, 0) - /C(mw, 0) + 2/C(m„„, 0) + 3e?^„,,,(0) + 3e?,„,,„(0) ] 

OTT 

+ ^AA^ [ 5/C(m„,, A) + /C(m„„, A) + /C(m,„, A) + /C(m,„, A) + 2/C(m,d, A) + 2/C(m,rf, A) 
+26;,,,,,(A) + 26;,,„,„(A) - 46;,„,,,(A) ] ) , (B4) 
where Grja,vbi^) = ^r?a%(^("^r,,, A), /C(m^,, A), /C(mx, A)), Hn^^n, is given by 



+ 7 — o ;r^7 — o ^ 



the mass, mx, is given by vn?^ = \ + ^fi)^ ^'^'^ ^("^tt, ^) is defined in eq. flT^ (where 
now the m^r dependence is made exphcit). Note that rriab refers to the Goldstone-boson mass 
with quark content a and b (hence m.„± = mud, etc.); j and / label the sea quark masses. 
The finite-size corrections to the neutron mass are given by 

+3e^,„*(o) + 66;,„,,,(o) + 36;,„,,(o) ] 

+ [ 2/C(mrfrf, 0) + /C(m,,, 0) + /C(m;rf, 0) - /C(m,rf, 0) + SQ^^^^M + ^S^u,vM ] 

OTT 

+ [ 5/C(m„rf, A) + }C{mM, A) + /C(m,d, A) + }C{mid, A) + 2/C(m,„, A) + 2/C(mi„, A) 
yvr 

+2^,,,,, (A) + 2^,„,,„(A) - 4^,„,,,(A) ] ) . (B6) 
In the isospin limit, one has 

2 2 

5lM^ = -^^[/C(m.,0) + 8/C(m;,0)] - ^[/C(m.,A) + /C(m:,A)] 

+ {591 + ^9 a) [ /C(m^, 0) - /C(m:, 0) ] , (B7) 

where we have used the fact that ^^^ ^^^(A) — |/C(m7r, A) in the isospin limit. Here 
denotes the mass of a pion made of one valence quark and one sea quark. These expressions 
further collapse down to isosp in-symmetric QCD in the limit — > m^. 



2. Nucleon Magnetic Moments 

The most general charge matrix whose matrix elements reduce to those of QCD (keeping 
the valence-quark charges fixed) is fiij 

Q^^'^) ^ diag [+\^ Qj , qi, q, , qi) ■ (B8) 
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The finite-size corrections to the proton magnetic moment in PQQCD are 



+ 



9\ 
9 



4 y{muu. 0) - 5 y{mud. 0) - 4 y{mju. 0) - 4 y{miu, 0) 



9 



y{mud, 0) + y{muu, 0) - y{mju, 0) - y{miu, 0) 



9l 



+ 1^ [ 3^(m„,, 0) + 4 y{m^^, 0) - 3 3^(m,,, 0) + 3 3^(m,d, 0) 



-4 3;(m,„, 0) + 3 3^(mirf, 0) - 4 y{miu, 0) 



+9i 



y{mju, 0) - 3^(m„„, 0)] + 



9i9a 



y{m,u, 0) - y{muu, 0)] 



6 



3 3 
y{mju, 0) - 3^(m„„, 0) + - y{mjd, 0) - - 3^(m„rf, 0) 



+qi 



''^91 



3^(mz„,0)- 3^(m„rf,0)] 



6 



3^(mz„, 0) - y{m^d, 0) + ^ 3^(m,d, 0) - ^ y[mM, 0) 



+ ^ [ Km,,, A) - 3^(m„„,A)-6y(m„,,A)- ^K,, A) 



3^(m,-„,A)- 3^(m,,,A)+ 3^(mi„,A) 



+ - g,- ( yira^u. A) + 2 3^(m„,, A) - 3^(m,„, A) - 2 3^(m,,, A)) 
+ ^ * ( 3^('^^.d, A) + 2 3^(mrfrf, A) - y{mi^, A) - 2 3^(mzrf, A)) 



,(B9) 



where 3^(m7r, A) is defined in eq. (j25|l (where now the m^r dependence is made explicit). 
The finite-size corrections to the neutron magnetic moment are 



+ 



9i9a 



9\ 
9 



7 3^(m„,, 0) + 2 3^(m;,, 0) + 2 3^(m,,, 0) - 2 3^(m,,, 0) 



9 



y{m^d. 0) + 3^(m;,, 0) - y{mud. 0) - 3^(m,,, 0) 



+ [ 3 3^(m„„, 0) + 2 y(m„rf, 0) - 3^(m,,, 0) + 3^(m,,, 0) 
io 



-3 3^(m,„, 0) + y{mid, 0) - 3 3^(m,„, 0) 



+9i 



y{m,d. 0) - 3^(m,rf, 0)] 



+ 



9\ 



y{mjd, 0) - y{m^d. 0) + ^ 3^(m,„ 0) - ^ 3^(m„„, 0) 



3^(mw,0)- 3^(m,,,0)] 



^7? 



3^K,, 0) - y{mdd, 0) + ^ 3^(mi„, 0) - ^ 3^(mw, 0) 



+ 



54 



y{mdd, A) - 4 3^(m„„, A) + 9 3^(m„,, A) - y{m,d, A) 



13 



+4 y{mju, A) - y{mid, A) + 4 y{miu, A) 
+3 Qj ( y{mud, A) + 2 3^(m,„, A) - y{mjd, A) - 2 3^(m,„, A)) 
+3 qi ( :^(mdd, A) + 2 :^(m„d, A) - y{mid, A) - 2 :^(m,„, A)) ] ) . (BIO) 

In the isospin limit (with qj = qi = 0), one has 

hfc - ^ (^^f[y(m^,0) + 8y(ml,0)] + ^^y(m^,A) 

- (51 + ^9a) [ y{m^: 0) - y{ml, 0) ]) % . (Bll) 

These expressions further collapse down to isospin-symmetric QCD in the limit — > ttIt^. 
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